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Abstract 

The onset of chaos in one-dimensional spinning particle models derived from pseudoclassical 
mechanical hamiltonians with a bosonic Duffing potential is examined. Using the Melnikov method, 
we indicate the presence of homoclinic entanglements in models with general potentials for the 
spins, and thus show that chaotic motions occur in these models. 
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Suzuki and Maeda [1] has recently showed that the motions of spinning particles in the 
Schwarzschild black hole spacetime could be chaotic. This is an important result due to its 
relevancy to the detection of gravitational waves from black hole coalescences , and it 
has also aroused our interest in considering the chaotic motions of spinning particles. 

On the other hand, spinning particle motions have been considered in 0, 0] from a 
different point of view. The authors there described spinning particles as particles moving 
in a supersymmetric space, which they called the spinning space, with Grassmann-valued 
vector variables in addition to the ordinary spacetime coordinates. Hence, the interplay 
between chaos and spinning particle models involving pseudoclassical mechanics jf| looks 
rather intriguing. 

In this letter we would like to consider the simplest one-dimensional spinning particle 
model (Eq. below) derived from pseudoclassical mechanical systems . This is just the 
spinning particle model considered by Berezin and Marinov Q], or the so-called G3 model 
according to the classification scheme in We hope that the investigation of this simple 
model will shed some light on the chaotic behaviors of more complicated spinning particle 
models in higher dimensions or even in curved spaces. 



to detect the 



Here we use an analytic method called the Melnikov technique 
onset of chaotic motions. Suppose in an integrable system there are homoclinic orbits em- 
anating and terminating at the same unstable fixed point. Time-dependent perturbations 
are then added to this system. In the perturbed system the stable and unstable orbits split. 
An integral called the Melnikov function evaluated along the unperturbed homoclinic orbit 
measures the transversal distance between the perturbed stable and unstable orbits on the 
Poincare section. The presence of isolated zeros in the Melnikov function indicates com- 
plicated entanglements between the two perturbed orbits and thus the presence of chaotic 
behaviors. 

We start with the bosonic part of our model. We choose the one-dimensional Duffing 
hamiltonian, 

H=£ + V(x), (1) 

where 

V(x) = -x 2 + x A . (2) 
The equations for the homoclinic orbit (separatrix) with energy E = are particularly 
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simple, 

x s (t) = sech(v / 2t), (3) 
Pa (t) = -V2 sech(v / 2t) tanh(v / 2t). (4) 

This is the main reason we choose the Duffing potential. On the other hand, we expect 

similar analysis as the one we carry out below can be applied to other potentials with 

homoclinic or heteroclinic orbits. 

To consider pseudoclassical models with Grassmann variables, we concentrate on those 

with three Grassmann variables £1, £2, and £3, 

2 • 3 

H = P - + V{x)- 1 - J2 (5) 

i,j,k=l 

where Wi(x) are the potentials for the Grassmann variables. The hamiltonian can be written 



as 
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H = ^- + V(x) - W(x) ■ S (6) 



1 

Si = --< ijkCjCb- (7) 



where the "spin" S is defined by 
The dynamical equations are 

dec / _ \ 

* = ^ (8) 

f = -V'(x) + W'(x) ■ S, (9) 
^ = -W(x) x S, (10) 

where the prime ' denotes d/dx. Note that from here on, we work directly with this set of 
dynamical equations without further reference to the Grassmann nature of S as defined in 
Eq. (0). That is, we do not restrict 

5? ~ 0. (11) 

Hence, we are dealing with the spinning particle models defined by the dynamical equations 
rather than the general pseudoclassical systems with Grassmann variables In effect we 
just use the pseudoclassical hamiltonians as a convenient and systematic means to derive 
the dynamical equations for the spinning particles. 
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FIG. 1: Time evolution of S± in Eqs. (|14|) - (|16j) with the initial conditions as given in Eq. (|17|). 
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FIG. 2: Time evolution of 52 in Eqs. (|14|) - (|16j) with the initial conditions as given in Eq. (|17|) . 

To detect the onset of chaotic motions in the spinning particle models above, we treat 
the spin S as a perturbation. That is, we assume that 



x,p~0(l) and Si ~ 0(e) 



(12) 



for some small parameter e. Then the additional spin terms in Eq. Q have the form of 
time-dependent perturbations which may trigger chaotic behaviors. 

To make the discussion more concrete we first take a simple form for the potentials W(x), 



W 1 = x; W 2 = l; W 3 = 0. 



(13) 



We shall return to the consideration of more general W(x) below. To obtain the lowest 
order solutions of Sj, we just replace W(x) in the spin dynamical equations (Eq. (JTUJ) ) by 
W(x s (t)), that is, 

dSi 



dt 



(14) 
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FIG. 3: Time evolution of S3 in Eqs. (|14|) - (|16j) with the initial conditions as given in Eq. (|17|). 



dS 2 
~dt 
dS 3 
dt 



sech(V2t)S 3 

Si - sech(V2t)S 2 



(15) 
(16) 



where we have written Si = eSi. We can numerically solve this set of differential equations 
quite easily using, for example, Mathematica. The solutions with the initial conditions 



Si(0) = l ; S 2 (0) = S 3 (0) = 0, 



(17) 



are plotted in Figs. [T]|21 Si and S3 are oscillatory, while S 2 is mostly constant except for a 
small region near t = 0. Note that the choice of the initial conditions here is quite arbitrary. 
We have checked the behaviors of the spins for other sets of initial conditions like Si(0) = 0; 
S 2 (0) = 1; 6*3(0) = and Si(0) = 0; ^(O) = 0; S 3 (0) = 1, and similar results like those in 
Figs. HH21 are obtained. 

It is well known that external oscillatory perturbations may induce chaotic behaviors. 
We can see that this is indeed the case here by calculating the Melnikov function 0,0], 



M(t 



dt Ps (t + t )W'(x s (t))- S(t) 

"diV2 sech[V2(t + to)] tanh[V2(t + t )]Si(t), 



(18) 



where to is used to parametrize the location along the unperturbed homoclinic orbit. M(to) 
is plotted in Fig. |3J The Melnikov function measures the transveral distance on the Poincare 
section between the perturbed stable and unstable orbits passing through the unstable fixed 
point. The infinite number of zeros as seen in Fig. 0] indicates the entanglement of the 
two orbits, and this signals the occurrence of chaotic behavior for the perturbed homoclinic 
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FIG. 4: Melnikov function of the spinning particle model in Eq. © with the potentials for the 
spins W(x) as given in Eq. (|TT?|) . 

orbits. When the perturbation terms, or the spin terms here, become larger and larger, the 
chaotic behavior will spread to other parts of the phase space. 

For more general W(x) potentials, we can argue that similar chaotic behaviors as those 
with the simple potential in Eq. (fT^j) discussed above will again occur. In this case the 
Melnikov function is still given by Eq. ()18)) but with general W(x). The key point to the 
argument is to consider the large to behavior of M(to). This is just the region near the 
unstable fixed point where complicated homoclinic entanglements may occur. Now in the 
integrand in Eq. (JTSJ), 

p s (t + t ) = -V2 sech[V2(t + *„)] tanhfv^t + t )]. (19) 

which is exponentially small except for 

t ~ -t - (20) 

Hence, to estimate M(t ) for large t , we need only to consider the behaviors of W'(x s (t)) 
and S(t) for t ~ — 1 or for large |t|. 

For large |t|, x s (t) is also exponentially small. Suppose, for small x, 

W i ~a i + b i x + ---, (21) 

where a« and bi are constants. Then, for large |t|, Wi can be approximated by a iy and the 
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spin equations become 



d 
dt 



s 2 

\S 3 J 



a 3 — a 2 
— a 3 ai 

y a 2 -oi o y 



(22) 



a? + a| 



ai The solutions 



The eigenvalues of the matrix above are and ±ia, where a = 
of Si are therefore linear combinations of sin(at), cos(at), and some constants. Except for 
the case when all — 0, some of the components Si must be oscillatory. 
Similarly, for large |t|, 

W'i ~ h. (23) 

Hence the Melnikov function for large to can De written approximately as a linear combina- 
tion of 



/oo 
dtp s {t + t ) = 0, (24) 

f°° iia ( ira \ 

M2{to) — dt p s (t + t ) sin (at) = -= sech — -= cos(at ), (25) 

J-oo V2 \2a/2/ 

/OO yj-Q / \ 

dt p s (t + to) cos(at) = — —j= sechf -j= J sin(at ), (26) 

where we have used the fact that p s (t) is odd in t. The Melnikov function M(t ) for general 
W(x) potentials with the small x behavior as shown in Eq. (|21j) is therefore an oscillatory 
function with infinity number of simple zeros at large to- This is sufficient to indicate the 
occurrence of chaotic behaviors for these general spinning particle models. 

In summary, we have indicated, using the method of Melnikov functions, the presence of 
homoclinic entanglements and thus the onset of chaos in one-dimensional spinning particle 
models with a Duffing potential and some general potentials for the spins. We expect similar 
results for other potentials with homoclinic or heteroclinic orbits. Moreover, the models we 
consider in this letter are derived from one-dimensional pseudoclassical systems with three 
Grassmann variables, so we also expect chaos to occur for spinning particle models derived 
from systems with three or more Grassmann variables where more degrees of freedom are 
involved. 

As we have mentioned at the beginning, chaotic behaviors do occur for the motions of 
spinning particles in a Schwarzschild black hole spacetime which can be described by a 



super symmetric theory. This shows that we should extend our investigation of chaos in 
spinning particle models to higher dimensions as well as to models in curved spacetimes. 
We plan to do that in our future publications. 
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the National Science Council of the Republic of China under contract number NSC89-2112- 
M-032-021. 
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